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Abstract 

Symmetry is an important feature of many con- 
straint programs. We show that any symmetry act- 
ing on a set of symmetry breaking constraints can 
be used to break symmetry. Different symmetries 
pick out different solutions in each symmetry class. 
We use these observations in two methods for elim- 
inating symmetry from a problem. These methods 
are designed to have many of the advantages of 
symmetry breaking methods that post static sym- 
metry breaking constraint without some of the dis- 
advantages. In particular, the two methods prune 
the search space using fast and efficient propa- 
gation of posted constraints, whilst reducing the 
conflict between symmetry breaking and branch- 
ing heuristics. Experimental results show that the 
two methods perform well on some standard bench- 
marks. 

1 Introduction 

Symmetry occurs in many problems. For instance, certain 
workers in a staff rostering problem might have the same 
skills and availability. If we have a valid schedule, we may 
be able to permute these workers and still have a valid sched- 
ule. We typically need to factor such symmetry out of the 
search space to be able to find solutions efficiently. One pop- 
ular way to deal with symmetry is to add constraints which 
eliminate symmetric solutions (see, for instance, 1 11 |2j [3j 
m [51 El |3)- A general method is to add constraints which 
limit search to the lexicographically least solution in each 
symmetry class. Such symmetry breaking is usually simple 
to implement and is often highly efficient and effective |8; 
|9 1 . Even for problems with many symmetries, a small number 
of symmetry breaking constraints can often eliminate much 
or all of the symmetry. 

One problem with posting symmetry breaking constraints 
is that they pick out particular solutions in each symmetry 
class, and branching heuristics may conflict with this choice. 
In this paper, we consider two methods for posting symmetry 
breaking constraints that tackle this conflict. The two meth- 
ods exploit the observations that any symmetry acting on a 
set of symmetry breaking constraints can be used to break 
symmetry, and that different symmetries pick out different 



solutions. The first method is model restarts which was pro- 
posed in |10|. We periodically restart search with a new 
model which contains a different symmetry of the symmetry 
breaking constraints. The second method posts a symmetry of 
the symmetry breaking constraint dynamically during search. 
The symmetry is incrementally chosen to be consistent with 
the branching heuristic. Our experimental results show that 
both methods are effective at reducing the conflict between 
branching heuristics and symmetry breaking. 

2 Background 

A constraint satisfaction problem (CSP) consists of a set of 
variables, each with a domain of values, and a set of con- 
straints specifying allowed combinations of values for sub- 
sets of variables. A solution is an assignment to the variables 
satisfying the constraints. We write sol{C) for the set of all 
solutions to the constraints C. A common method to find a 
solution of a CSP is backtracking search. Constraint solvers 
typically prune the backtracking search space by enforcing a 
local consistency property like domain consistency. A con- 
straint is domain consistent iff for each variable, every value 
in its domain can be extended to an assignment satisfying the 
constraint. We make a constraint domain consistent by prun- 
ing values for variables which cannot be in any satisfying as- 
signment. During the search for a solution, a constraint can 
become entailed. A constraint is entailed when any assign- 
ment of values from the respective domains satisfies the con- 
straint. For instance, Xi < Xn is entailed iff the largest value 
in the domain of Xi is smaller than the smallest value in the 
domain of Xn- A constraint is dis-entailed when its negation 
is entailed. For instance, X <Y is dis-entailed if and only if 
the smallest value in the domain of X is larger than or equal 
to the largest value in the domain of Y. 

CSPs can contain symmetry. We consider two common 
types of symmetry (see 111 II for more discussion). A variable 
symmetry is a permutation of the variables that preserves so- 
lutions. Formally, a variable symmetry is a bijection cr on the 
indices of variables such that if Xi = di, . . . , Xn — dn is a 
solution then ^cr(i) — di, . . . , = dn is also. A value 

symmetry, on the other hand, is a permutation of the values 
that preserves solutions. Formally, a value symmetry is a bi- 
jection 9 on the values such that if Xi = di, . . . , Xn — dn is 
a solution then Xi — 9{di), . . . , Xn = d{dn) is also. Sym- 
metries can more generally act on both variables and values. 



Our methods also work with such symmetries. As the inverse 
of a symmetry and the identity mapping are symmetries, the 
set of symmetries of a problem forms a group under compo- 
sition. We will use a simple running example which has a 
small number of variable and value symmetries. This exam- 
ple demonstrates that we can use symmetry itself to pick out 
different solutions in each symmetry class. 

Running Example. The all interval series problem (probOO? 
in CSPLib.org hl2^ ) asks for a permutation of to n — 1 
so that neighbouring differences form a permutation of 1 to 
n — 1. We model this as a CSP with Xi = j iff the ith number 
is j, and auxiliary variables for the neighbouring differences. 
One solution for n = II is: 

Xi,X2,...,Xn = 3,7,4,6,5,0,10,1,9,2,8 (a) 

This model has a number of different symmetries. First, there 
is a variable symmetry Urev that reverses any solution: 

Xi,X2,...,Xn = 8,2,9,1,10,0,5,6,4,7,3 (b) 

Second, there is a value symmetry Oinv that inverts values. If 
we subtract all values in ^ from 10, we generate a second 
(but symmetric) solution: 

Xi,X2,...,Xn = 7,3,6,4,5,10,0,9,1,8,2 (c) 

Third, we can do both. By reversing and inverting (toj), we 
generate a fourth (but symmetric) solution: 

Xi,X2,...,Xn = 2,8,1,9,0,10,5,4,6,3,7 (d) 

The model thus has four symmetries in total: aid ( the identity 

mapping), arev, Oinv, and Oinv ° Orev * 

3 Symmetry breaking 

One common way to deal with symmetry is to add con- 
straints to eliminate symmetric solutions 111. Two important 
properties of symmetry breaking constraints are soundness 
and completeness. A set of symmetry breaking constraint 
is sound iff it leaves at least one solution in each symmetry 
class, and complete iff it leaves exactly one solution. 

Running Example. Consider again the all interval series 
problem. To eliminate the reversal symmetry a^ev, we can 
post the constraint: 

< Xn (1) 

This eliminates solution © as it is the reversal of ([op. To 
eliminate the value symmetry Oinv, we can post: 

Xi<5, Xi = 5 ^ ^2 < 5 (2) 

This eliminates solution ((cp as it is the inversion of ((op. Fi- 
nally, to eliminate the third symmetry Oinv ° (^rev where we 
both reverse and invert the solution, we can post: 

{Xi,...,Xe) <ie. (10-Xn,...,10-X6) (3) 

This eliminates solution ((op as it is the reversal and inversion 
of (O. Note that of the four symmetric solutions given ear- 
lier, only Q with Xi — 2 and Xn = 7 satisfies all these 
symmetry breaking constraints. The other three solutions are 
eliminated. Thus (|7} fo (O are a sound and complete set of 
symmetry breaking constraints. Jit 



We now show that any symmetry acting on a set of sym- 
metry breaking constraints itself breaks the symmetry in a 
problem. Different symmetries pick out different solutions 
in each symmetry class. To prove this, we need to consider 
the action of a symmetry on a symmetry breaking constraint. 
Symmetry has been defined acting on assignments. We lift 
this definition to constraints. The action of a variable sym- 
metry on a constraint changes the variables on which the con- 
straint acts. More precisely, a variable symmetry a applied 
to the constraint C{Xj, . . . , Xk) gives C{Xcr(j) , ■ ■ ■ , Xcr(fe)). 
The action of a value symmetry is also easy to compute. A 
value symmetry 9 applied to the constraint C{Xj, . . . ,Xk) 
gives C{0{Xj),...,e{Xk)). 

Running Example. To illustrate how we can break symmetry 
with the symmetry of a set of symmetry breaking constraints, 
we consider symmetries of (Q, ((21) and 

If we apply arev to (Q, we get an ordering constraint that 
again breaks the reversal symmetry: 

This simplifies to: 

Xn < Xi 

If we apply arev to (O, we get constraints that again breaks 
the inversion symmetry: 

Xn < 5, Xn = 5 ^ Xw < 5 

Finally, if we apply arev to Q, we get a constraint that again 
breaks the combined reversal and inversion symmetry: 

{Xn,...,Xe) <ie. (10-Xi,...,10-X6) 

Note that of the four symmetric solutions given earlier, only 
((cp satisfies arev o/dZ]), (H)) and 0. 

We can also break symmetry with any other symmetry of 
the symmetry breaking constraints. For instance, if we apply 
Oinv ° o'rev to (O, we get fl Constraint that again breaks the 
reversal symmetry: 

10 - Xn <10-Xi 

This simplifies to: 

Xi < Xn 

If we apply Oinv ° f reu to (O, we get constraints that again 
breaks the inversion symmetry: 

10-Xii<5, 10-Xii = 5 ^ 10-Xio < 5 

This simplifies to: 

Xn > 5, = 5 ^ Xio > 5 

Finally, if we apply Oinv ° o'rev to we get a constraint that 
again breaks the combined reversal and inversion symmetry: 

(10-Xn,...,10-X6) <,ex {Xi,...,Xe) 

Note that of the four symmetric solutions given earlier, only 
((op satisfies Oinv ° Orev o/dZ]), (O and (O. ^ 

The running example illustrates that we can break sym- 
metry with a symmetry of a set of symmetry breaking con- 
straints. We now prove that this holds in general: 



Any symmetry acting on a set of symmetry breaking 
constraints itself breaks symmetry. 

More precisely, if a set of symmetry breaking constraints is 
sound, then any symmetry of these constraints is also sound. 
Similarly, if a set of symmetry breaking constraints is com- 
plete, then any symmetry of these constraints is also com- 
plete. 

Observation 1. Given a set of symmetries S ofC, if S is a 
sound (complete) set of symmetry breaking constraints for E 
then cr(S') for any a £ Y, is also a sound (complete) set of 
symmetry breaking constraints for S. 

Proof: (Soundness) Consider s G sol{C U S). Then 
s G sol{C) and s S sol{S). Hence cr(s) G sol{C). 
Since s G sol{S), it follows that a{s) G sol{a{S)). Thus, 
cr(s) G sol{C U cr(5)). Hence, there is at least one solution, 
(t(s) in every symmetry class of C U cr{S). That is, cr{S) is a 
sound set of symmetry breaking constraints for S. 

(Completeness) Consider s G sol{C U cr{S)). By a similar 
argument to soundness, a^^{s) G sol{C U S). Hence, there 
is at most one solution in every symmetry class of C U cf{S). 
That is, cr{S) is a complete set of symmetry breaking con- 
straints for S. □ 

Different symmetries of the symmetry breaking constraints 
pick out different solutions in each symmetry class. Thus, 
if the branching heuristic is going towards a particular so- 
lution, there is a symmetry of the symmetry breaking con- 
straints which does not conflict with this. 

Observation 2. Given a symmetry group E, a sound set S 
of symmetry breaking constraints for E, and any complete 
assignment A, then there exists a symmetry a in E such that 
A satisfies a{S). 

Proof: Since the set of symmetry breaking constraints S is 
sound, it leaves at least one solution (call it B) in the same 
symmetry class as A. That is, B satisfies S. Since A and B 
are in the same symmetry class, there exists a symmetry cr 
in E with (j{A) = _B. E forms a group so also contains the 
inverse symmetry a^^. Since B satisfies S and <j{A) = B, 
it follows that a{A) satisfies S. Hence a^^{(T{A)) satisfies 
a-^{S). That is, A satisfies a^^iS). □ 

4 Model restarts 

We start with a simple application of these observations. 
To tackle conflict between branching heuristics and symme- 
try breaking constraints. Heller et al. propose using model 
restarts IIIOI . In this method, backtracking search is restarted 
periodically, using a new model which contains different 
symmetry breaking constraints. By posting different sym- 
metry breaking constraints, we hope at some point for the 
branching heuristic and symmetry breaking not to conflict. 
Our observations that any symmetry acting on a set of sym- 
metry breaking constraints can be used to break symmetry, 
and that different symmetries pick out different solutions, 
provide us with precisely the tools we need to perform model 
restarts to any domain (and not just to interchangeable vari- 
ables and values as in IIIOI ). When we restart search, we 
simply post a different symmetry of the symmetry breaking 
constraints. We experimented with several possibilities. The 



simplest is to choose a symmetry at random from the sym- 
metry group. We also tried various heuristics like using the 
symmetry most consistent or most inconsistent with previous 
choices of the branching heuristic. However, we observed the 
best performance of model restarts with a random choice of 
symmetry so we only report results here with such a choice. 

Running Example. Consider again the all interval series 
problem and posting symmetries of the symmetry breaking 
constraints (O, (|2]l and The following table gives the 
amount of search needed to find an all interval series of size 
n — 11 using a branching heuristic that branches in order 
on the variables introduced to represent neighbouring differ- 
ences in the series, trying values in numerical order. This 
would seem to be a good branching heuristic since it can find 
a solution without backtracking. 



Symmetry posted 
of O to © 


Branches 


Time to solve/s 




1 


0.00 




222,758 


13.74 




425,765 


24.99 




170,425 


10.23 



It is clear from this table that the different symmetries of 
the symmetry breaking constraints interact differently with 
the branching heuristic. In particular, the identity symme- 
try does not confiict in any way as the branching heuristic 
goes directly to the following solution at the end of the first 
branch: 

Xi,X2,...,Xii = 0,10,1,9,2,8,3,7,4,6,5 

This solution is consistent with aid o/dl}. (0 <^nd (O. 

The other symmetries of the symmetry breaking constraint 
confiict with the branching heuristic. In particular, the fol- 
lowing symmetry breaking constraints confiict with this so- 
lution: <7rev o/dTJ as Xn = 5 > Xi = 0, 9inv o/ ([7} as 
10 — Xi = 10 > 10 — Xii = 5 and Oinv ° (^rev o/0 as 
10 — Xii = 5 > Xi = 0. As a result, posting these sym- 
metries of the symmetry breaking constraints increases the 
search needed to find a solution. 

Model restarts will help overcome this confiict. Suppose 
we restart search every 100 branches and choose to post a 
random symmetry of (Q, (|2) and Let t be the average 
number of branches to find a solution. There is j chance that 
the first restart will post aid o/dTJ, (|2]l and (O. In this situa- 
tion, we find a solution after 1 branch. Otherwise we post one 
of the other symmetries of (O, (|2| and We then explore 
100 branches, reach the cutoff and fail to find a solution. As 
each restart is independent, we restart and explore on aver- 
age another t more branches. Hence: 

t = ii + ^(ioo + i) 
4 4 

Solving for t gives t = 301. Thus, using model restarts with 
a random symmetry of the symmetry breaking constraints, we 
take just 301 branches on average to find an all interval series 
of size n = 11. Jit 

Note that posting random symmetries of the symmetry 
breaking constraints is not equivalent to fixing the symmetry 



breaking and randomly branching. As we saw in the example, 
different symmetries of the symmetry breaking constraints 
interact in different ways with the problem constraints. Al- 
though the problem constraints are themselves initially sym- 
metrical, branching decisions quickly break the symmetries. 

5 Posting constraints dynamically 

We now consider a more sophisticated use of the observa- 
tions that any symmetry acting on a set of symmetry break- 
ing constraints itself breaks symmetry, and that different sym- 
metries pick out different solutions in each symmetry class. 
We will incrementally and dynamically post a symmetry of 
the symmetry breaking constraints which is consistent with 
the branching decisions made so far. Thus, if the branching 
heuristic is smart or lucky enough to branch immediately to a 
solution, symmetry breaking will not interfere with this. 

Running Example. Consider again the all interval series 
problem. Suppose we begin by trying Xi = 10. Since the Xi 
are all different, Xn S [0, 9]. Hence, the symmetry breaking 
constraint Xn < Xi is entailed. This is a rev o/dZJ. It is 
also 9inv of (O}. We do not yet need to commit to which of 
these two symmetries of the symmetry breaking constraints 
we will post. We are sure, however, that we are not posting 
Gid or Oinv ° cfrev of the Constraints (|7} to di]). These two 
symmetries would require Xi > Xn, and this is dis-entailed. 
We therefore post Xn < Xi and continue search. Jit 

In the example, we posted symmetry breaking constraint 
once they are entailed. When there are only a few symme- 
tries, we can easily implement this with non-backtrackable 
variables and reification. Suppose we reify the two ordering 
constraints: 

Bi ^ (Xi < Xn), B2 ^ (Xn < X^) 

We then make the Boolean variables, Bi and B2 non- 
backtrackable so that, once they are instantiated, their value 
remains on backtracking. We assume that our solver posts the 
conclusion of an implication when its hypothesis is entailed. 
Suppose Xi < Xii is entailed. Then Bi will be set true. As 
Bi is non-backtrackable, Xi < Xu will be posted. Unfor- 
tunately, posting symmetry breaking constraints like this as 
soon as they are entailed may be a little eager. 

Running Example. Consider again the all interval series 
problem. As before, suppose the branching heuristic has set 
Xi = 10, and we have posted the entailed symmetry break- 
ing constraint Xu < Xi. Now Xi > 5 is also entailed. This 
is Oinv of the first inequality in (|2]l. If we post this, we com- 
mit to breaking symmetry with Oinv of (|7} to (O. However, 
this would rule out breaking symmetry with o^ev o/(I2} 
which are also still consistent with the branching decisions so 

far- 
Suppose we next branch on X\\ = 5. The assignments to 
Xi and Xii are only consistent with Oinv o/dU and o/dij- 
In fact, both of these constraints are now entailed. However, 
Xi = 10 and Xu = 5 are not consistent with posting ffrev 
o/dU- This would require that: 

{Xii,...,Xe) <iex (10-Xi,...,10-X6) 



This is dis-entailed. Hence, our branching decisions have 
committed us to break symmetry with Oinv of d2} to dU- We 
therefore post these constraints. If search continues, we will 
discover the unique solution consistent with symmetry break- 
ing and the initial branching decisions: 

Xi,X2,...,Xii = 10,0,9,1,8,2,7,3,6,4,5 

There is a tradeoff between posting symmetry breaking 
constraints early (so propagation prunes the search space) and 
late (so we do not conflict with future branching decisions). 
We propose the following rule for when to post symmetry 
breaking constraints. The rule only posts symmetry break- 
ing constraints once the branching heuristic has forced their 
choice. It would, however, be interesting to explore other 
more eager or lazy rules. Suppose S* is a set of symmetry 
breaking constraints for E, and we have posted T, a symme- 
try of a subset of S. A symmetry cr G E is consistent with 
T iff T is entailed by cr(5) and inconsistent otherwise. A 
symmetry cr e E is eliminated by posting some symmetry 
breaking constraint c iff a is consistent with T but inconsis- 
tent with T U {c}. The forced symmetry rule is defined as 
follows: 

Given a set of symmetry breaking constraints, if 
during backtracking search a symmetry of one of 
these constraints is entailed, this symmetry is con- 
sistent with previously posted symmetry breaking 
constraints, and all symmetries eliminated by this 
entailed constraint are inconsistent with the current 
state then we post the entailed constraint. 

We first show that this rule is sound. 

Observation 3. Given a set of symmetries S of C, if S is 
a sound set of symmetry breaking constraints for E then the 
forced symmetry rule using S is a sound symmetry breaking 
method. 

Proof: The rule only permits constraints of a particular sym- 
metry to be posted. By Observation 1, this is sound. □ 

In general, this rule may not be complete even when given 
a complete set of symmetry breaking constraints. However, it 
is easy to modify the rule so that it is complete. Whenever we 
reach a solution, we simply pick a consistent symmetry and 
post all the symmetry breaking constraints associated with 
this symmetry. We can also define a common property of 
many symmetry breaking constraints for which the unmodi- 
fied rule is complete. A set of symmetry breaking constraints 
S for the symmetries E of C is proper iff S is sound and 
complete for E and every non-identity symmetry in E maps 
any solution of S* U C onto a different solution. With a proper 
set of symmetry breaking constraints, each solution within a 
symmetry class is associated with a different symmetry. For 
instance, constraints dTJ to d3]l form a proper set of symmetry 

We new prove that with a proper set of symmetry breaking 
constraints, the forced symmetry rule is a sound and complete 
symmetry breaking method. That is, it will find exactly one 
solution in each symmetry class. 

Observation 4. Given a set of symmetries E of C, if S is a 

proper set of symmetry breaking constraints for E then the 
forced symmetry rule is both sound and complete. 



Proof: (Soundness) Immediate as a proper set is sound. 

(Completeness) Consider the first solution visited. As the 
set of symmetry breaking constraints is proper, only one sym- 
metry of these constraints will be entailed. All other symme- 
tries are inconsistent with the current state and are eliminated. 
The forced symmetry rule therefore post this symmetry of 
the symmetry breaking constraints. By Observation 1, as the 
symmetry breaking constraints are complete, this eliminates 
all other solutions in the same symmetry class. □ 

Finally, we observe that with certain symmetry breaking 
constraints, the forced symmetry rule is equivalent to posting 
symmetry breaking constraints as soon as they are entailed. 
For symmetry breaking constraints like Xi < Xn, as soon 
as the constraint or its negation is entailed, all variable sym- 
metries are either consistent or they are eliminated. 

6 Interchangeable variables and values 

To test these two symmetry breaking methods, we consider 
a common type of symmetry where variables and values par- 
tition into interchangeable sets 1 13; 14 1. This is sometimes 
called piecewise variable and value symmetry. We chose this 
class of symmetry other the many other types of symmetry 
studied in the past as it was used in the previous experimen- 
tal study of model restarts II 1011 . Suppose that the n variables 
partition into a disjoint sets and variables within each set are 
interchangeable. Similarly, suppose that the m values parti- 
tion into h disjoint sets and values within each set are inter- 
changeable. We will order variable indices so that ^p(i) to 
Xp(i+i)-i is the ith partition of variables for \ < i < a, and 
value indices so that d^Q) to dq(j+i)-i is the jth partition of 
values for 1 < j <h. 

Flener et al. II 141! proved that we can eliminate the expo- 
nential number of symmetries due to such interchangeability 
with a polynomial number of symmetry breaking constraints: 

^p(i) < • • • < -^p(j+l)-l 

Gcc([Xp(j), . . . , Xp(j+i)_i], [di, . . . , dm], [Oi, . . . , O^J) 

• ■ ■ 'O^O)) >lox ■ ■ • >lcx (Oj(j + i)_i, • ■ ■ + 

Where i e [1, a] and j G [1,6], and Gcc counts the number 
of occurrences of the values in each equivalence class of vari- 
ables. That is, O] = \{k\Xk = dj,p{i) < k < p{i + 1)}|. 

The signature of dk is (O^, . . . , O^), the number of occur- 
rences of dk in each varialDle partition. The signature is in- 
variant to the permutation of variables within each equiva- 
lence class. By ordering variables within each equivalence 
class, we prevent permutation of interchangeable variables. 
Similarly, by ordering the signatures, we prevent permutation 
of interchangeable values. 

We will post symmetries of these symmetry breaking con- 
straints. We consider symmetries that act along two degrees 
of freedom: the order of interchangeable variables within a 
variable partition, and the order of the signatures of inter- 
changeable values within a value partition. Let a be some 
permutation of the indices of interchangeable variables. Then 
we can break the symmetry of variable interchangeability 
with the following symmetry of the variable ordering con- 
straints: 



Similarly let 6 be some permutation of the indices of inter- 
changeable values. Then we can break the symmetry of value 
interchangeability with this symmetry of the signature order- 
ing constraints: 

7 Experiments 

We used model restarts and the forced symmetry rule to 
post symmetries dynamically of the symmetry breaking con- 
straints of Flener et al. 1 14|. Problems are coded into Gecode 
2.2.0. We evaluated the two methods on the same two bench- 
mark domains used in previous studies of symmetry break- 
ing for interchangeable variables and values |15|. Experi- 
ments were run on an 2-way Intel Xeon with 6MB of cache 
and 4 cores in each processor running at 2GHz. All in- 
stances were terminated after 10 minutes. We used smallest 
domain as a variable ordering heuristic in each experiment. 
For value ordering heuristic, we used lexicographical, anti- 
lexicographical and random orderings. 

Our experiments are designed to test two hypotheses. The 
first hypothesis is that these two methods are less prone to 
conflict between branching heuristics and symmetry break- 
ing. The second hypothesis is that these two methods explore 
a smaller search tree than dynamic methods like SBDS. This 
is due to both the propagation of the posted symmetry break- 
ing constraints and the need to limit SBDS to just generators 
to make it computationally tractable. We limit our compari- 
son of dynamic methods to comparison against SBDS. Whilst 
there is a specialized dynamic symmetry breaking method 
for interchangeable variables and values, experiments in 1 10| 
show that this is several orders of magnitude slower than 
static methods. In addition, dominance detection methods 
like SBDD are shown to be three orders of magnitude slower 
than static methods in 1 10|. Finally, we used SBDS to break 
just generators of the symmetry group as breaking the full 
symmetry group quickly ran out of memory. We used SBDS 
with two different sets of generators: one set has a genera- 
tor that exchanges each pair of consecutive variables/values 
in each partition; the other has a generator that exchanges the 
first two variables/values in each partition and one that rotate 
all variables/values. We got similar results with both and thus 
here report only results for the first set, denoted SBDS-pair in 
Table 1. 

The first set of experiments uses random graph coloring 
problems generated in the same way as the previous exper- 
imental study in I15II . There is a variable for each vertex 
and not-equals constraints between variables corresponding 
to connected vertices. All values in this model are inter- 
changeable. In addition, we introduce variable symmetry by 
partitioning variables into interchangeable sets of size at most 
8. We randomly connect the vertices within each partition 
with either a complete graph or an empty graph, and choose 
each option with equal probability. Similarly, between any 
two partitions there is equal probability that the partitions 
are completely connected or independent. Results for graphs 
with 40 vertices are shown in the top half of Table 1 . 

The second set of experiments uses a more structured 
benchmark which is again taken from a previous experimen- 
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Table 1: Static vs Dynamic posting of symmetry breaking 
indicates optimality was not proved), "t" is the runtime in 



constraints on Graph Coloring and Concert Hall Scheduling, "opt" is the quality of the solution found (* 
seconds, "b" is the number of backtracks. The best method for a problem instance is in bold font. 



tal study II15I . In the concert hall scheduling problem, we 
have n applications to use one of m identical concert halls. 
Each application has a start and end time as well as an offer 
for the hall. We accept applications so that their intervals do 
not overlap and the profit (the sum of the offers of accepted 
applications) is maximized. We randomly generate instances 
so that applications are split into partitions of size at most 8 
and within each partition all applications have the same start 
and end time and offer Our model assigns Xi = j if the i*^ 
application is accepted and placed in hall j, and Xi ^ m + 1 
if it is rejected. Variables corresponding to applications in the 
same partition are interchangeable. Values divide into two 
partitions: the values 1 to ni are interchangeable, while the 
value m + 1 is in a separate partition. Results for instances 
with 40 applications and 10 halls are shown in the bottom half 
of Table 1. 

The results support both our hypotheses. Both methods 
are less prone to conflict between symmetry breaking and the 
branching heuristic. With both SBDS and our forced sym- 
metry rule for dynamically posting symmetry breaking con- 
straints, we obtained the same results with the lexicographical 
and the (inverse) anti-lexicographical value ordering heuris- 
tic. With model restarts, results with the lexicographical 
and the anti-lexicographical value ordering heuristic are suf- 
ficiently similar that we only report the former. Our second 
hypothesis, that the two methods explore a smaller search tree 
than SBDS is also confirmed. SBDS was unable to prove op- 
timality in all but one instance. In addition, on the harder 
graph coloring instances, both methods tend to outperform 
the static method. It is hard, however, to choose between 
the two methods. The model restarts method offers slightly 
better performance on the graph coloring instances, whilst 
our method of dynamically posting static symmetry break- 
ing constraints offers better performance on the concert hall 
scheduling instances. 

8 Other related work 

Closest in spirit to our forced symmetry rule for dynami- 
cally posting symmetry breaking constraints is SBDS |16l 
[l2l[T8 |. SBDS can work with any type of branching decision 
but for simplicity we assume that branching decisions are of 
the form Var = vol. All current implementations of SBDS 
make this assumption. If we have a symmetry cr, the partial 
assignment A and have explored and rejected Var = vol 
then on backtracking, SBDS posts: 

ct{A) a{Var ^ val) 

This ensures that we never explore the symmetric state to 
the one that has just been excluded. Our forced symmetry 
rule also posts static symmetry breaking dynamically dur- 
ing search. However, the two methods differ along three 
important dimensions. First, SBDS posts symmetry break- 
ing constraints when backtracking and exploring the second 
branch of the search tree; the forced symmetry rule, on the 
other hand, can post symmetry breaking constraints down ei- 
ther branch. Second, SBDS posts symmetries of the current 
nogood; the forced symmetry rule, on the other hand, can 
post any type of symmetry breaking constraint. Here, for in- 
stance, it posts ordering constraints on the signatures. Third, 



whilst neither method conflicts with the branching heuristic if 
the branching heuristic goes directly to a solution, the forced 
symmetry rule may conflict with the branching heuristic later 
in search. Constraint propagation on constraints posted by the 
forced symmetry rule can prune values that branching might 
have taken. 

There are a number of other related methods. Jefferson 
et al, have proposed GAPLex, a hybrid method that also 
combines together static and dynamic symmetry breaking 
|:19|. However, GAPLex is limited to dynamically posting 
lexicographical ordering constraints, and to searching with 
a fixed variable ordering. As a consequence, GAPLex per- 
forms poorly when there are large numbers of symmetries. In 
addition, GAPLex is unable to profit from effective dynamic 
variable ordering heuristics. Puget has also proposed "Dy- 
namic Lex", a hybrid method that dynamically posts static 
symmetry breaking constraints during search which works 
with dynamic variable ordering heuristics [20] . This method 
adds lexicographical ordering symmetry breaking constraints 
dynamically during search that are compatible with the cur- 
rent partial assignment. In this way, the first solution found 
during tree search is not pruned by symmetry breaking. Un- 
fortunately Dynamic Lex needs to compute the stabilizers of 
the current partial assignment. This requires a potentially 
expensive graph isomorphism problem to be solved at each 
node of the search tree. Whilst Dynamic Lex works with dy- 
namic variable ordering heuristics, it assumes that values are 
tried in order Finally Dynamic Lex is limited to posting lex- 
icographical ordering constraints. This is problematic when 
there are many symmetries. A direct comparison of our meth- 
ods with Dynamic Lex would be interesting but poses some 
challenges. For instance. Heller et al. II 1011 do not compare 
model restarts with Dynamic Lex, arguing: 

"It is not clear how this method [Dynamic Lex] can 
be generalized, though, and for the case of piece- 
wise variable and value symmetry, no method with 
similar properties is known yet." 

9 Conclusions 

We proved that any symmetry acting on a set of symmetry 
breaking constraints itself breaks symmetry, and that differ- 
ent symmetries pick out different solutions in each symmetry 
class. These observations can be used to reduce the conflict 
between symmetry breaking and branching heuristics. We 
have studied two methods for breaking symmetry that tackle 
this conflict. The first method uses model restarts which was 
proposed in |10|. We periodically restart search with a new 
model which contains a random symmetry of the symmetry 
breaking constraints. The second method posts a symmetry of 
the symmetry breaking constraint dynamically during search. 
The symmetry is incrementally chosen to be consistent with 
the branching heuristic. The two methods benefit from prop- 
agation of the posted symmetry breaking constraints, whilst 
reducing the conflict between symmetry breaking and branch- 
ing heuristics. Experimental results demonstrated that the two 
methods perform weU on some standard benchmarks. 
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